arXiv:1506.09130vl [astro-ph.HE] 30Jun2015 


The Pulsating Pulsar Magnetosphere 

K.H. Tsui 

Instituto de Fisica - Universidade Federal Fluminense 
Campus da Praia Vermelha, Av. General Milton Tavares de Souza s/n 
Gragoatd, 24-210-346, Niteroi, Rio de Janeiro, Brasil. 

tsui@if .uff .br 

ABSTRACT 


Following the basic p rinciples of a charge separated pulsar magnetosphere 
(iGoldreich fc Juharull96911 . we consider the magnetosphere be stationary in space, 
instead of corotating, and the electric held be uploaded from the potential distri¬ 
bution on the pulsar surface, set up by the unipolar induction. Consequently, the 
plasma of the magnetosphere undergoes guiding center drifts of the gyro motion 
due to the transverse forces to the magnetic held. These forces are the electric 
force, magnetic gradient force, and held line curvature force. Since these plasma 
velocities are of drift nature, there is no need to introduce an emf along the held 
lines, which would contradict the = E ■ B = 0 plasma condition. Furthermore, 
there is also no need to introduce the critical held line separating the electron 
and ion open held lines. We present a self-consistent description where the mag¬ 
netosphere is described in terms of electric and magnetic helds and also in terms 
of plasma velocities. The helds and velocities are then connected through the 
space charge densities self-consistently. We solve the pulsar equation analytically 
for the helds and construct the standard steady state pulsar magnetosphere. By 
considering the unipolar induction inside the pulsar and the magnetosphere out¬ 
side the pulsar as one coupled system, and under the condition that the unipolar 
pumping rate exceeds the Poynting hux in the open held lines, plasma pressure 
can build up in the magnetosphere, in particular in the closed region. This could 
cause a periodic openning up of the closed region, leading to a pulsating magne¬ 
tosphere, which could be an alternative for pulsar beacons. The closed region can 
also be openned periodically by the build-up of toroidal magnetic held through 
a positive feedback cycle. 


Subject headings: (stars:) pulsars : general 
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1. Rotating Magnetosphere 


The basic model of a pulsar had long been recognized as a perfect c onductor (n eutron 
star) rotating with an angular velocity f2o in a hiper magnetic held B flGoldl 1196811 . Con¬ 
sidering the rotational axis and the magnetic axis both aligned in the same sense, charged 
particles in the neutron star are driven by the magnetic force to the surface according to their 
signs, with electrons towards the polar regions and ions towrds the equatorial region. This 
is the Faraday unipolar induction that operates homopolar generator, which extracts the en¬ 
ergy of a conducting rigid rotator. Due to its rotation in the presence of a magnetic held, the 
conducting pulsar surface is not an equipotential surface, contrary to the static case. If the 
pulsar were surrounded by a vacuum, the separated charges would accumulate on the pulsar 
surface by their signs, and a counter electric held would build up inside the neutron star to 
counteract the magnetic force. A steady state surface charge distribution or potential distri¬ 
bution would be reached when the two forces inside the pulsar cancel out, and there would 
be a quadrupole electric held in space surrounding the neut ron star. The unipo l ar ind uction 
would cease to operate then. However, the seminal work of iGoldreich fc Julian! (119691) noted 
that such a vacuum conhguration was unstable, and there had to be a pulsar magnetosphere 
to channel the surface charges. As a result, the counter electric held inside the pulsar would 
be small, and the unipolar induction would continue to operate. The steady state potential 
distribution on the pulsar surface would be uploaded to the equipotential magnetic held lines 
of the magnetosphere. Quite often, held line breaking near the pulsar is invoked to provide 
an electron-positron background pair plasma to hll the magnetosphere. 


The equation of motion of a charge particle in the magnetosphere is given by 


P 


dv 

dt 


= pqE + J X B — Vp — upv + pg 


With the plasma inertia on the left side and plasma pressure, collision, gravitational force on 
the right side neglected, comparing to the hiper magnetic and electric helds, the macroscopic 
structure of the magnetosphere is determined on ly by the balance between electric and 


magn etic forces, which is the force-free held model flScharlemann fc Wagonerll 19731: 
1974 ) with 


Okamoto 
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pqE + J X 5 = 0 . 


( 1 ) 


Although plasma inertia and plasma pressure are neglected in the macroscopic description, 
these terms in the equation of motion give rise charged particle gyro motion about magnetic 


field lines and transverse dri 


t motions acros s field lines in the microscopic single particle 


Thomoson 

1962; 

Schmidt 

1966; 

Jackson 

1975) 


197511 ■ As a result of these drifts, 


macroscopic drift currents and plasma flows, corresponding to the guiding center motions, 
are generated in Eq. [H as will be addressed in Sec.3. 


With the surface charges having ready access to the field lines, plasma density will be 
high enough in the magnetosphere such that plasma conductivity a tends to be infinite, the 
Ohm law reads 


E + v X B = —J = 0 . (2) 

cr 

For a quasi-neutral normal magnetohydrodynamic (MHD) plasma, the plasma velocity v in 
Eq. [2] amounts to the center of mass velocity. For a charge separated plasma, this equation 
coincides with Eq. [1] other than an overall fator of the space charge density pq. Nevertheless, 
this does not mean that Eq. [2] is redundant, since it describes the magnetosphere in terms 
of velocities while Eq. [T] does it in terms of fields. These two descriptions are then coupled 
through the space charge densities. 

Probably induced by the lighthouse scenario of pulsar beacons, the pulsar magnetosphere 
is often interpreted as such that the magnetic held lines are rigidly anchored on the pulsar 
surface. The entire magnetosphere thus corotates with the pulsar at its angular velocity, and 
the charge separated plasma is being dragged along with it. Although for an aligned rotator 
and under axisymmetry, it is indistinguishable whether or not the magnetosphere is rotating. 
Nevertheless, in an oblique rotator, it could generate the lighthouse beacons and magnetic 
dipole radiation. In this scenario, the rotating magnetic held B drags the plasma generating 
a rigid rotor velocity v. Applying this rotating plasma velocity v back to a stationary frame 
in space, and considering the magnetic held B be stationary, which is indistinguishable from 
a rotating one, Eq. [2] then gives the magnetospheric electric held E in the stationary frame. 
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We remark that only under the axisymmetric aligned rotator case that the magnetic held 
can be considered stationary and/or rotating by option. So the plasma velocity v is an 
input parameter in Eq. [21 the driving energy is the rotating magnetic held, and the electric 
held E is the response which is called the rotation induced electric held. Some authors 
introduce a non-rotation induced longitudinal held — VE by taking into consideration of the 

he light cylinde r 


interst ellar plasma potentiah or by considering the centrifugal outhow at 
(LC) f Mestel et ah 1985 : 


Fitzpatrick fc Mestellll988 : 


Goodwin et ah 


29941 : 


Timokhin 


299d). 


Other authors have considered a very stringent supply of positive charges from the pulsar 
surface, due to probably the surface physics of neutron star, and proposed that vacuum gaps 
could be formed separating the charg ed regions from the pulsar sur f ace and/or separating 
the c harged regions from each other (IRuderman fc Sutherlandlll975 : 


Jackson 


1976 


Michel 


1979h . 


In this corotating magnetic held model, the unipolar induction in the neutron star is 
clearly decoupled from the magnetosphere. The magnetic dipole radiation, the angular mo¬ 
mentum outhow of the stellar winds, and the current loop in the presence of —VE are 
attributed to the braking torque on the pulsar indirectly through the general energy con¬ 
servation principle. However, the use of the energy conservation principle implies a direct 
coupling with the neutron star, which is not rehected in this scenario. As a matter of fact, if 
the magnetosphere held lines are anchored on the pulsar surface, the (dipole) magnetic held 
inside the pulsar will be rotating as well. Consequently, there is no relative rotation between 
the pulsar and the magnetic held inside. For this reason, there will be no unipolar induction 
as well. 


As a matter of fact, a steady stat e magnetosphe re decoupled from the pulsar implies that 
the emission problem of the beacons fjMelroselll978l) can be constructed from the equilibrium 
magnetosphere. Although the validity of separating equilibrium and emission into two diher- 
ent issues cannot be disproved, this approach could have over simplihed the pulsar dynamics. 
Recently, a time dependent description of magnetosphere where the plasma equations are 
coupled to the Maxwell equations to account for induction ehects. This broader approach 


has been adopted to study the formation of m agnetosphere current sheets ([Komissarov 


2996 


Spitkovsky 


2996; 


Tchekhovskoy et al.[ [2913h. and the line ar instabilities in magnetosphere 


that could be responsable for emissions fjUrpin|[2912[. [29141) . 
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2. Neutron Star-Stationary Magnetosphere Coupled System 


To overcome the rotation ambiguities of plasma velocity and magnetic field, we consider 
the magnetic field B be stationary in space, and the electric field E is uploaded to the 
equipotential held lines of the magnetosphere according to the potential distribution on the 
pulsar surface estblished by the unipolar induction. This means that the charge separated 
magnetospheric plasma has a negligible contribution to the electric held there. In this case 
of stationary magnetic held B, the driving energy is the uploaded electric held E which is 
the input parameter in Eq. [2l and the plasma velocity v is the response. The magnetosphere 
is directly coupled to the unipolar induction of the pulsar. The charge separation how 
inside the pulsar generates a current. The interaction of this current with the magnetic held 
produces a direct braking torque in addition to the contributions of dipole radiation and 
stellar winds which slow down the rotation rate. Finally, we remark that magnetic held is a 
real physical quantity, while magnetic held line is only a physical concept derived from the 
held line equation to help visualize the action of magnetic held. It is rather unphysical to 
consider held lines being rigidly attached to the pulsar surface corotating with it. The notion 
of a rotating magnetosphere is probably cultivated in the oblique rotator case, where the 
magnetosphere is seen to rotate obliquely about the pulsar rotational axis. As obliqueness 
goes to zero, we then have an aligned rotator. For an oblique rotator, we choose to phrase 
the magnetosphere as wobbling (not rotating) about the rotational axis. Consequently, as 
the obliqueness goes to zero, the magnetosphere stops to wobble and becomes stationary in 
space. 


Here, in this paper, we reconsider the iGoldreich fc JulianI fjl969l ) analysis of charge 
separated plasmas. In particular, there are two inconsistencies in their work. The hrst in¬ 
consistency is the space charge density expression which is incompatible with the curvature 
of the closed ion magnetic field lines near the pulsar. To resolve this inconsistency, they 
resorted to a corotating electron cloud above the closed magnetic held lines. The second one 
is the introduction of the far zone quasi-neutral interstellar plasma potential to the near zone 
pulsar magnetosphere to draw plasma outhows in the open held lines. The held line that has 
the potential of the interstellar plasma is the critical held line with no charge how. This in¬ 
terstellar emf is inconsistent to the E\\ = 0 plasma condition in the magnetosphere. Actually, 
many publications in recent decades have concentrated on the pulsar equation itself without 
mentioning neither the question of space charge density nor the condition of the critical held 
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line. These simplifications are eqnivalent to considering the pulsar magnetosphere as filled 
with a quasi-neutral normal plasma instead of a charge separated plasma. 

The charge separated plasma in this stationary magnetosphere undergoes an E x B 
drift which happens to have an angular velocity of the pulsar, instead of being dragged along 
by the corotating magnetosphere. As for the plasma outflow, we consider it be driven by 
plasma drifts, not by interstellar emf. Finally, with a stationary magnetosphere in space, 
coupled to the unipolar induction in the pulsar interior, we offer an alternative pulsar beacon 
mechanism by suggesting a pulsating magnetosphere. To understand this dynamic picture, 
we assume that the unipolar induction is delivering more power to the magnetosphere than 
the open field lines can drain. This leads to an energy built-up in the closed field lines 
represented by plasma pressure. Although plasma pressure is neglected in Eq. [H it appears 
indirectly through plasma drifts, which could open up the closed field lines near LC. 

Observationally, radio pulsars are known for their stable repetitive pulse profiles. How¬ 
ever, the pulse profile that characterizes each radio pulsar is obtained as the average over 
many many shots. From shot to shot, the profiles are far from identical. Therefore, the 
period of oblique rotator is not well defined from shot to shot, which is not supposed to be 
in the lighthouse scenario. On the other hand, such shot to shot variations are expected in 
a pulsating magnetosphere. 


We review the basic concepts of charge separated pulsar magnetosphere and the plasma 
velocity v of Eq. |2] in Sec.3. The space charge density pq is derived in the presence of a 
poloidal and toroidal field, and in the presence of drift velocities in Sec.4. Under the con¬ 


cept of guiding center drifts, t 


le idea of a critical field line is reconsidered in Sec.5. We 


then 




bl 


ow the earlier works of 


StnrrockI (1197lh. IScharlemann fc Wagoner! (119731 ). 


Beskin fc Malvshkin 


(bd), 


Goodwin et ah 


Okamoto 


(120041) . to solve the magnetosphere of 
an isolated pulsar analytically via the force-free pulsar equation. Nevertheless, we take a 
different approach to split the pulsar equation into two coupled equations through a sep¬ 
arating function. Each one of the coupled equations is solved by separation of variables 
for the global structure of the pulsar magnetosphere in Sec.6. The standard monopole-like 
steady state magnetosphere is presented in Sec.7. A self-consistent description of magneto¬ 
sphere in terms of fields, velocities, and space charge densities is outlined in Sec.8. Under 
continuous unipolar pumping, the standard magnetosphere is generalized to the pulsating 
magnetosphere in Sec.9, and with conclusions in Sec. 10. 
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3. Plasma Velocity 

It is important to remark that Eq. |5] is written in an inertial frame where the plasma 
moves with velocity v with respect to the stationary magnetic and electric helds. In the 
plasma moving frame, the Lorentz transformed electric held is E' = E + v x B and the Ohm 
law reads E' = 0 there. This Ohm law in the limit of inhnite a is known as the frozen held 
condition in MHD plasmas, and happens to be identical to Eq. [H and we have E • n = 0 
and E\\ = E ■ B = As a matter of fact, for a charge separated plasma, Eq. [1] and Eq. [2] 
amount to the same equation always. However, Eq. [1] describes the magnetosphere in terms 
of helds, while Eq. [2] considers the magnetosphere in terms of plasma velocities. The helds 
and the velocities are then connected through the space charge density. 

Here we call to the attention that the condition of E\\ = 0 allows an arbitrary longitudinal 
velocity U|| along the magnetic held line. Consequently, the plasma velocity can be projected 
in components parallel and perpendicular to the magnetic held as 


U = U|| + Uj_ , 

= ve = ^(^j- X B) . 


The frozen held plasma condition shows that ue is given by the E x B guiding center drift of 


the gyro motion du e to the transverse part of the electric held E flThompson 


1962 


Schmidt 


19661: lJacksonlll975l) . under the constrain of < c which implies the electric held strength 
be E± < cB for a magnetized plasma. This ve drift is charge independent, and for a 
quasi-neutral normal plasma, this drift generates no current and represents only the huid 
drift velocity. But for charge separated plasmas, the E x B drift corresponds a current how. 
Consequently, we have the picture of a stationary pulsar magnetosphere with plasma drifting 
with respect to the magnetic held. 

In cylindrical coordinates {r,(j),z), the axisymmetric magnetic held can be represented 


as 


f Id^ 1 , ld^\ 

( — ) ’ 

\ r oz r r or J 


B — Aq^VT X V0 “1“ q:AV0) — Aq 


(3) 












IIqJ — Aq 


dA T 2 dA 

-a—,-V^^H-,+a 


dz 


r dr 


dr 


where A^ carries the dimension of the poloidal magnetic flux such that H/ is a dimensionless 
flux function, and a = A^jA q represents the relative ratio of the toroidal flux to the poloidal 
flux. The function A is related to the axial current through 


27raAoA = po4 • (4) 

From the magnetic held line equation, it is easy to show that the poloidal held lines are given 
by the contours of T 


^(r,z) = C . (5) 

The fact that the longitudinal electric held vanishes means that any potential diherence 
along the magnetic held line will be shorted out by plasma mobility. With the longitudinal 
electric held E\\ = 0 and the magnetic held B represented by Eq. [3l the transverse electric 
held in the magnetosphere can be expressed through the equipotential held lines of Eq. [5] as 


Ei_ = -EoVT , 

where Vq is the equator-pole voltage drop on the pulsar surface, and T is assigned by the 
corresponding surface potential label. 

Alternatively, we can also project the plasma velocity in poloidal and toroidal compo¬ 
nents by writing v = Vp + and B = Bp + B^. From Eq. O we get 0 = —= Vp x Bp by 
axisymmetry, and —Ep = Vp x B^p + x Bp. We therefore have respectively 



Vp KyBp , 
( ^v^(i} T ^0) ^ • 


( 6 ) 

( 7 ) 
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Loading the potential distribution of the pulsar surface onto the equipotential magnetosphere 
poloidal field lines, we have 


— Ep = rVL^cf) X Bp , 

where r measures the radial position of a point r in the magnetosphere. Substituting this 
result to Eq. [7] then gives and therefore 


V = Vp + V(j, = k^B + • (8) 

We call to the attention that drifts with an angular velocity of the pulsar. This equation 
is often interpreted as the ’beads on a wire’ model, where the plasma (beads) slides with a 
velocity KyB^f, along a rotating field line (wire). We have rederived this equation with 
a stationary magnetosphere where the rotation comes from the guiding center drift. 

We have expressed the plasma velocity in two orthogonal projections. The first one is 
with respect to the magnetic field lines, and the second one is with respect to the cylindrical 
coordinate system. These two projections are independent to each other. Only when the 
magnetic field is purely poloidal with B^ = 0 that the two expressions are correspondent, 
and we can identify 


nil Vp K^Bp , 

= W = he = rQo(f) . 
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4. Space Charge Density 


In the magnetosphere, the electric field E is not the only transverse force that causes 
plasma drift. Magnetic field line curvature and field gradient are other two major transverse 


Schmidt 


1966 


forces. Therefore, besides the vp drift, we a. 
f Thompson 1962 : 


Jackson 


so have the curvature and gradient B drifts 


19751) given by 


VRc = 

VyB = 


2e|| 1 Rc 0. 


q 52 


VB 


xB) , 


(9a) 

(9b) 


where e\\ = muy/2 and e± = mv‘]_/2 are the charge particle energies parallel and perpendicular 
to the magnetic field. Therefore, these drifts are mass dependent as well as field dependent. 
In a thermal plasma, these energies are measued by the plasma pressure. Although plasma 
inertia and pressure are neglected in Eq. [H they appear indirectly through these two drift 
velocities. The direction of these two drifts depends on the sign of charge q. Even in a 
quasi-neutral normal plasma, these are current generating drifts, while the electric drift ve 
represents the plasma fluid velocity only. In a charge separated plasma, all three drifts 
generate currents. These plasma drifts are driven by transverse forces, and they are not 
emf driven velocities. To understand this, we recall that for a uniform magnetic field, charge 
particles follow a circular gyro motion along the magnetic field lines according to the equation 
of motion due to the plasma inertia. In the presence of a transverse force to the magnetic 
field such as an electric field, a centrifugal force due to curvature, and an inhomogeneous 

be t r ansformed int o a displacing cycloid 
1^: I Jackson! Il975l) , such that the gyro 


magnetic field, etc., t he circular gyro motion wi 
across the field lines flThompsonl Il962 ; 


Schmidt 


guiding center presents a transverse drift. Although Eq. [2] gives the guiding center electric 
drift directly, the gyro motion itself is due to the plasma inertia. Therefore, Eq. [2] actually 
includes the plasma inertia implicitly. We remark that these drift velocities can be very 
large and become relativistic. Should the electric field E be time varying, there would be a 
polarization drift as well which we will ignore here. 


We should point out that in Eq. [2] only the ve drift is represented on the left side. The 
other two current generating drifts, that are contained in J on the right side of Eq. [21 do not 
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appear in Eq. [2] because the plasma conductivity a is infinite. For this reason, the plasma 
drift velocity is given by 


V drift — {vrc + VVb) + Ve ■ 

Since the plasma drift velocities are obtained from the equation of motion of each charged 
plasma, the plasma drift velocity vdrift is beyond the scope of Eq. |2l as discussed in Sec.l, in 
the sense that the MHD expression of n = v\\+ve differs from vdrift because n|| ^ {vrc+v\/b)- 

The space charge density, which connects the fields of Eq. [1] and the velocities of Eq. |2l 
in the magnetosphere is, by Eq. [2l 


Pq = -eoV • {ve X B) = eo{vE ■ PoJ - B - V x ve) ■ 

Writing J = PqVdrift = Pqi^E + vrc + v\ib)-i considering a general poloidal and toroidal fields, 
and representing ve by the projection of Eq. [H]to evaluate the —B - V xve term on the right 
side, we get 


[Pq - eoPoivE - k^B) • J] = [1 - eo/no(TE - t^vB) ■ Vdrift]Pq = -2eo5 • flo • (10) 

In the particular case of a purely poloidal field, we can identify ve = rVL^cj), and Eq. [TOl 
becomes 


{Pq - eo/^ohe • J) = (1 - eo/iohe • Vdrift)pq = -2eo5 ■ file 


Goldreich fc Julian! fjlQGQI) had derived a space charge density (their Eq.8) with J = pqVE 
alone that reads 


(1 — eofLQV%)pq — —2eoB ■ Gq • 
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In the absence of vrc and hys, the left side of the above equation looks like the relativistic 
Lorentz factor, which is just coincidental. This equation predicts <Q for closed ion held 
lines, which is consistent only to the outer portion of the dipole-like closed ion held lines. 
For the inner portion with Bz > 0, they postulated a negative charged cloud corotating with 
the held lines. 

We antecipate the result in Sec.8 that < 0, such that the {ve — t^vB) term is positive. 
Near the pulsar, vrc and v^b are large, where the plasma pressure is very high, on the 
left side and we could have Bz > 0. Farther out from the pulsar with plasma pressure 
decreasing we have Bz < 0. Since vrc and Fys drifts are mass dependent, they would 
be dominant again as the closed held lines approach LC, due to their relativistic inertia. 
We would have > 0 on approaching LC. For this reason, the closed region of a steady 
state magnetosphere can only extend to some intermediate radius, staying away from LC. 
Although plasma pressure and plasma inertia are not considered in Eq. [1] explicitly, plasma 
inertia appears implicitly through the guiding center drifts of the gyro motion and plasma 
pressure through the coefficients of the gradient B and curvature drifts. We have given a 
qualitative understanding of these ehects in the closed region through the drift velocities, 
and we will return to this issue in Sec.8. 
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5. Critical Field Line 

In terms of cylindrical coordinates, the toroidal components of Eq. [H which describes 
the magnetosphere in terms of fields, gives 


{aVA X V0) X (VIE' x V0) = 0 , 

A = A(d/) , (11) 

whereas the poloidal components yield the pulsar equation 


1 - i—Y —— + a^A{m) ; = 0 

tl J r or aW 


( 12 ) 


where = c/VLq is the LC radius. Because of the singular coefficient of the highest 
derivative which vanishe s at r/ rj, = 1, Eq. [12] is a singular equation. As indicated by 


Scharlemann fc Wagoned (119731 ) in Sec.Ill, the steady state magnetosphere, which is sta¬ 


tionary in space, is solved within LC, and is then analytically continued beyond LC where 
the pulsar equation is not valid with the tengential velocity = rhlo > c. Since plasma 
parameters like plasma pressure and velocities are not considered in the pulsar equation, 
analytic continuation of the field structure beyond LC amounts to solving the pulsar equa- 


tion in that region with s o 


Contoponlos et ah! (119991 ). 


utions on both s ides matched across LC, as done numerically by 


TimokhinI (1200(11 ). and others. 


Analytic solutions of thi s equation had been extensive 


were not obtained explicitly (jScharlemann fc Wagoner! 1 1973 : 


y studied, yet g lobal solutions 


Okamoto 


I 974 I) . So far, solu¬ 


tions are developed by numerical iterratiq n s, under the boundary condition of a dipole field a t 
the center no t ably by Contoponlos et ah ( 1999 ). Qgura fc Koiima ( 2003 ). Gruzinoy (12005 ). 


ContoDonlosI (120051 ) where the open and closed regions rotate differentially, and 


J2OO6I) 


Timokhin 


where the braking index is evaluated in terms of the magnetosphere dynamics. 


In the open field lines, electrons and ions are assummed to stream out along their 
respective open field lines with equal flux. These currents will be connected to form a return 
loop at a distant load (interstellar plasma). The electron and ion outflows are thought to be 
drawn by the action of the far zone quasi-neutral normal interstellar plasma floating potential 
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Einter = EVinter^^i wheie Vinter IS the potential difference between the far zone interstellar 
plasma and the near zone electron or ion open field lines, and $ is a dimensionless scalar 
potential fnnction. In particnlar, there will be a field line having the same potential of the 
interstellar plasma. This is the critical held line Tc, where there is no charge and no poloidal 
cnrrent how, that separates the electron and ion held lines. This means that dA{^)/d'^ = 0, 
and A(T) should be stationary at Tc. This critical held line should pierce LC with = 0. 

Nevertheless, using the interstellar emf to draw current hows in the open held lines 
is clearly in confrontation with the E\\ = 0 condition. This inconsistency can be removed 
should we note that the presence of a toroidal and a poloidal magnetic held in the open 
region can drive a poloidal and a toroidal ve drift respectively 


ve = Vp + v^ = -^Ee X + Bp) . 

The same can be applied to vjic and hys drifts. All these drift velocities will, therefore, 
produce a component along the magnetic held lines plus a toroidal rotation, as in Eq. [HI 
and there is no need to introduce neither the interstellar potential to the near zone nor the 
concept of a critical held line. 
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6. Analytic Pulsar Solution 

By considering the magnetic field B be stationary in space, and the electric field E be 
nploaded to the field lines from the pnlsar snrface, onr model has the plasma velocity v as 
the response in Eq. |2l With = rVL^cf), our model simply implies that the electric potential 
on the pulsar surface cannot be fully uploaded to the magnetosphere to all distances, even 
beyond LC. This means that the plasma conductivity a is finite in Eq. [2] on approaching LC, 
and collisional effect has to be taken into account, or equivalently, inertia effect. Likewise is 
in Eq. [H This should be the subject of further investigations. For the moment, we persuit 
Eq. [12] as it is. 

To construct the standard magnetosphere, we denote ^ = r/ri as the normalized radial 
coordinate and take a different approach to recast the pulsar equation as 




(1-a 


1 2d^ 

irl ^ 




M(T)' 

9T 


fm , 


(13) 


where we have required both sides of the equality be equal to a separating function /(T). 
We attempt to solve this equation analytically by choosing 


/(T) = , 


which renders 


2d^ 


'W 

(arL)M(T) 


M(T) 

(9T 


+ ’ 


(14) 


(15a) 

(15b) 


where = zjrL is the normalized axial coordinate. We write T(,^,<^) = R{^)Z{<,) by separa¬ 
tion of variables, and Eq. I15al gives 


d^Z 

'W 


{{krif — m?)Z = 0 , 


(16a) 
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d^R 


m^R , 


(16b) 


where m? is the separation constant. We have denoted the two constants in Eq. [H] and 
Eq. Il6bl by k'^ and m? respectively in quadratic form by choice. However, and m? are not 
necessarily be positive. As for Eq. Il5b( we take a linear A(\[') and get 


A(d') = -k,^ , 

9 rlT^ 

= [{cyKr,f + {kr,f{l-e)]R ■ 


(17a) 

(17b) 


Here, k^ has the dimension of an inverse scale length. We have chosen a minus sign explicitly 
in Eq. I17al to represent an overall downward axial current in the northern polar region. 
With the understanding that electrons and ions flow along their held lines by plasma drifts, 
A(\I') does not have to have a stationary point on LC which ables the linea r choi ce in 
Ea.ll7al Such choice of A(\ l/) = —kz^ was also used bv IScharlemann fc Wagoner! fjl973l ) and 
Beskin fc Malvshkin fjl998 h Our analytic approach to the pulsar equation has introduced 


a separating function /(d') = a separating constant m^, and an axial current function 
A(d') = -kz^. 


The Z{R) function can be solved readily to render 


^a(?) = Cae-^^ , (18a) 

Z,(<^) = , (18b) 


with = {{kri)"^ — m?) = (a| — m?) > 0, and 


= C'^cosn’q , (19a) 

Zl{q) = Oft sin kV , (19b) 
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with = —(m^ — a|) < 0. To solve for R{^), we make use of Eg. llTbI to cast Eg. I16bl 

as 


2-^ + [{{ak^rLf + {krif - 2w?) - {krLf^‘^]R 

= + [{al + al- 2m^) - ali^]R 

rP R 

= 2^^ + K - = 0 • (20) 

This eguation can be solved through a power series. Because of the boundary condition of 
z) = 0 on the polar axis r = 0, we should have R{^) = as ^ goes to zero. We 
therefore have 


^n. — 


R{0 = 

1 

0-2 — 


24 

1 


2 


^ 


■2(n + 2)(n+l)*d“”-2 


( 21 ) 


( 22 ) 


We note that Oq and Oi are the two independent constants of the second order differential 
eguation, Eg. |20l Other coefficients are generated by the recurssion formula, Eg. J22l The 
coefficients oq and oi generate even and odd power terms of the series respectively. Thus oq 
and oi are the multipliers of the even and odd series that we can write 


Ri.^y ®l) (^oReveni^: 1 ) 0 ) T ^iRoddi,^) 0 ) 1 ) • 


(23) 
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7. Standard Magnetosphere 

We note that the radial function R{^) is governed by the parameters (a^, a|) and the 
choice of the two independent coefficients (oq, Oi). With positive and a|, the coefficients 
ttn are alternating in sign, and R{^) could have oscillating solutions for small This func¬ 
tional form is suitable to construct the open ion held lines and the closed region of the pulsar 
magnetosphere. As ^ increases, R{^) becomes monotonically increasing, which is suitble for 
the split monopole structure far from the pulsar. As for the axial function Za{<^), it is a de¬ 
creasing function of <^. When this is combined with the radial function, together they could 
generate contours adequate for open held lines. For a given a^, we could scan a| over a given 
range to generate a set of held line morphologies. As a result, we represent the electron open 
held lines by taking the poloidal hux function as 


T(e,<^) = . (24) 

For simplicity, we consider the even terms of the power series by taking Oq = 1 and ai = 0. 
We choose to set = 20, say, for the radial solution. Taking a| = 8, the function R{^) is 
a monotonically increasing function, as is shown in Fig.l which is appropriate for electron 
open held lines. We should emphasize that the monotonic functional form of R{^) can be 
obtained under a wide range of parameters, not restricted to the chosen ones. Furthermore, 
we set = 5, for example, and we have the self-consistent parameters = 18, = 3. 

With = 5, the axial function Za{i;) is shown in Fig.2, and a set of poloidal held lines is 
shown in Fig.3. 

According to Eq. [3], the poloidal and toroidal magnetic 


Bp = Ao —X 0 

rLi 

^ . kzl ^ . kz I ^ 

= -aAo —-T = -A^ —-T 

Bp _ Ao 1 VT 

B^ ~ ^ 


helds are given by 

(25a) 

, (25b) 

(25c) 
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where the last equation gives the ratio of the two helds in the magnetosphere. However, this 
ratio is expressed in terms of their respective fluxes, Aq and A^p, without which a specific 
ratio cannot be determined. 

For the ion open held lines, we use the same representation of Eq. [2H and keep = 20 
and = 5. However, we change a| = 6 and other parameters become = 16, m? = 1. 
By changing a| = 8 to a| = 6, the radial function begins to present oscillations for small 
which is appropriate for the ion held lines. The functions R{C} and are given in Fig.4 

and Fig.5 respectively. We note that R{C} rises from zero and passes through a maximum 
and a minimum locally before it takes oh monotonically. The local maximum and minimum 
can become more pronounced by taking a smaller a\. A set of poloidal held lines is shown 
in Fig.6 which rehects the monopole-like solution of a charge separated magnetosphere in 
the near zone. In the boundary zone, the plasma tends to be quasi-neutral as it proceeds to 
meet the interstellar plasma in the far zone. The ^ < 1 part of the open held lines is the 
solution of the pulsar equation, and the > 1 part is the analytic continuation across LC. 

To construct the closed held lines of the ion plasma, known as the dead zone, we note 
that, with ion current outhow and electron return current outhow (closed at a distant load) 
along the open held lines equal in magnitude. Bp, = 0 in the closed region which requires 
al = 0. To represent the closed held lines, we need a radial function bounded by a root 
R{0 — 0 within LC. We also need a bounded axial function. We therefore write 


= m)Z',W) = C . (26) 

In order to have R{^) bounded by a root within LC, we choose oq = -l-l and oi = —1.5 
with diherent signs. Under the condition = 0, we take = —17, a| = 7, and = 5 
= —5), we have = 12, and the functions R{^) and are shown in Fig.7 and 

Fig.8 respectively. Once more, we recall that these parameters are neither unique ones nor 
restrictive. There is a good range in parameter space that reproduces the closed region. The 
boundary condition of = 0 on the equatorial plane at <^ = 0 is satished by A 

corresponding set of poloidal held lines is shown in Fig.9 with a diherent ^ scale of Fig.3 and 
Fig.6 for more clarity. In this hgure, the closed region terminates at an intermediate value 
of .^ = 0.80 and the reason of which has been discussed in Sec.4. The same hgure is redrawn 
with the same scale of Fig.3 and Fig.6, so that the three hgures can be superimposed on 
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each other to get the global features of the magnetosphere, as in Fig.10. This set of global 
held lines is a continuous function of the contour value C. More held lines could be added 
by assigning the in t ermed iate contour values C. Comparing to the numerical results, such as 


Contopoulos et ahl fjl999l) . our magnetospheric conhguration reproduces well the main global 


features. 
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8. A Self-Consistent Description of Magnetosphere 

So far, we have introduced plasma drift velocities to remove the inconsistency of the 
space charge density with the closed held lines, and the incinsistency of the E\\ = 0 condition 
with the emf driven plasma outhows. Furthermore, we have solved analytically the pulsar 
equation for the global structure of the magnetosphere. Above all, we have taken the view 
of a stationary magnetosphere in space under which the plasma drifts with a velocity v 
according to Eq. [2l As we have already mentioned, the magnetosphere is described in terms 
of helds by Eq. [H and in terms of velocities by Eq. [2l We should note that the force-free 
approximation is only valid for the held-description. As for the velocity-description, plasma 
inertia and plasma pressure are the prime factors. These two descriptions are then connected 
by Eq. [10] which describes the magnetosphere in terms of space charge densities. All these 
three aspects have to be self-consistent among each other. With a general poloidal and 
toroidal magnetic helds, the plasma velocities are given by Eq. jS] derived from plasma drifts. 
As for the space charge density, it is rehected by the B ■ Qq on the right side of Eq. [101 

Before addressing the consistency of the space charge density, let us hrst establish the 
magnitudes of the vrc and Fvs drifts by noting their ratio 


VRc _ ^ 
VVB Rc 


(27) 


where L is the scale length dehned by 1/L = S/B/B. Next, to compare with ve, we rewrite 
hh in terms of T which reads 


X B) = -l/o^(VT X B) . (28) 

For a poloidal magnetic held with T = C representing the held line, VT is directed inward 
towards the equatorial plane, and ve corresponds to a toroidal rotation of ve = rf2o0- For 
the Fvs drift of Eq. [Ob] we can evaluate Vi? by using the T = (F representation to get 


VB = Aoa^VT , 


(29) 
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where a has the dimension of an inverse scale length. Equivalently, has the dimension 
of a reference poloidal mangnetic field. As a result, Eq. [9b] reads 


X B) . (30) 

We note that this equation is identical in form for a positive charge q to Eq. |28l As for 
the vrc drift of Eq. jHa] we note that the radius of curvature Rc directs from the center of 
curvature outward to the held lines. Therefore, Rc and — VT are pointing outward, although 
they are not parallel. Consequently, both uys and vrc drifts are in the toroidal direction, as 
is ve- 

As for their amplitudes, based on Eq. [2H1 Eq. [301 and Eq. [OH all three velocities have 
the same 1/i?^ scaling, so their strengths can be compared by other factors. For vec, we 
note that Rc is small when it is near the pulsar and when it is near the apex on returning 
back to the conjugate point of the southern equatorial region. In between, the held lines are 
relatively hat, as shown in Fig.9, and Rc is large. As for v^b, the ratio of this velocity to ve 
is 


Ve _ qVo B 

A 2 * 

vvB Aoa^ 

Assumming the reference poloidal held comparable to B, these two velocities become 
comparable when the plasma energy (pressure) e± = qVo- The same can also be said for ey 
for vrc- This high plasma pressure can be found close to the pulsar surface, and it decreases 
as r increases. 

For a toroidal magnetic held in +0 direction, the ve drift will be poloidal but in the 
counter direction of closed held lines. Likewise are hye and vrc- By Eq. El we therefore have 
Ky < 0, and the factor {ve — KvB) in Eq. [TO] is positive. Should we consider the toroidal 
magnetic held be in —cj) direction, we would get Ky > 0, and the factor {ve — KyB) would 
be positive again. We remark that the space charge density of Eq. [TO] is valid for the entire 
magnetosphere, not just for the closed region. In the presence of poloidal and toroidal helds, 
the scalar product in the square bracket on the right side of Eq. [TO] has a poloidal and a 
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toroidal conribution. This space charge density has to be calculated together with the helds 
and plasma drifts in numerical solutions. 
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9. Pulsating Magnetosphere and Magnetosphere Instability 

We have constructed the standard magnetsphere in Sec. 7 where the closed region with 
= 0 is described by Eq. |26l On the other hand, with the self-consistent description of the 
magnetosphere based on fields, velocities, and space charge densities in Sec.8, there should 
be a in the closed region. On this question of a toroidal held, we note that, in a charge 
separated magnetosphere, a toroidal magnetic held can arise spontaneously. If we begin 
with a pure poloidal magnetic held in the polar region open electron held lines, equilibrium 
condition requires the electron plasma be at rest. This is an unstable conhguration, since any 
longitudinal movement of the electron plasma corresponds to a poloidal current that would 
generate a toroidal magnetic held. This toroidal magnetic held would enhance the poloidal 
drift velocity (current) itself. Since the held lines are open, an equilibrium conhguration of 
helds and currents could be reached. With pg < 0 for electron held lines, and with small 
due to small it is reasonable that the square bracket be positive and thus > 0. As 
for the ion open held lines, the toroidal held here begins to decline due to the ion outhows. 
With Pq > 0 for ion held lines, it is likely that we could have > 0 in some part of the held 
line and < 0 in other part. 

As for the ion held lines in the closed region, there still is a toroidal held due to the 
imperfect shielding of the electron axial current by the ion outhow in the open region. With 
Pq > 0, and with the presence of a toroidal held and a high plasma pressure near the 
pulsar, Eq. [10] could have B^ > 0 near the pulsar. Farther out with large the toroidal 
held becomes small, likewise is the plasma pressure, we could expect B^ < 0. With the 
poloidal held lines described by Eq. [26l the corresponding toroidal held in the closed region 
is represented by Eq. I25bl Because of the mirror symmetry between the upper (north) and 
lower (south) hemispheres, the toroidal helds have opposite sign in the two parts, which 
drives an outward radial current sheet on the equatorial plane. On reaching the apex of 
the closed region, this current sheet divides into two parts returning from the northern and 
southern closed poloidal held lines. We therefore have two current loops, one north and one 
south, to enhance the corresponding toroidal helds. This in turn enhances the radial current 
sheet and the current loops. In a closed hnite region, this positive feedback cycle can only 
be stablized in the presence of dissipations. Naturally, there is the usual toroidal current 
sheet on the equatorial plane beyond LC for the split monopole poloidal held lines. 

To hnd an alternative beacon mechanism, we hrst racall that the current model pictures 
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a corotating magnetic field that drags the magnetospheric plasma. In response to this, a 
rotation induced electric field is generated in the stationary frame of space. In this picture, 
the rotating magnetic field is the prime driver, and the magnetosphere is decoupled from the 
unipolar induction of the pulsar. Here, we view the unipolar induction inside the pulsar and 
the magnetosphere outside the pulsar as one coupled system. The unipolar induction delivers 
charges separated by their signs to the polar and equatorial regions at a given rate. The 
current flow inside the neutron star interacts with the magnetic field, generating a braking 
torque. The charges are then channelled to the magnetosphere along the field lines, and are 
ejected to the interstellar space as stellar winds through open held lines, which allows the 
unipolar induction to operate continuously. Considering the unipolar energy pumping rate 
exceeds the Poynting outhow rate along the open held lines, plasma density and pressure 
will build up in the magnetosphere, particularly in the closed region. In the outer part of 
the closed region extending to intermediate ^ with Bz < 0, hvs and vrc are small. However, 
as the plasma pressure build-up reaches there, these velocities begin to increase. Eventually, 
they get to the point of having > 0, hipping open the held lines and releasing the trapped 
plasma and magnetic energies. Furthermore, the held lines can also be hipped open in the 
outer part of the closed region by the positive feedback on the toroidal held which appears 
in the {ve — HvB) factor of Eq. [101 

Once depleted of plasmas, in the sense that plasma density is so low that Eq. [2] is not 
warranted, the magnetosphere will then be recharged all over to the standard conhguration, 
and the whole cycle starts anew. Instead of a stationary magnetosphere, we therefore have 
a scenario of a pulsating magnetosphere where the openning of the closed held lines either 
by unipolar pumping on drift velocities or by positive feedback on the toroidal held works 
as a magnetic switch, which ohers an altenative mechanism for pulsar beacons. In this 
alternative picture, an observer aligned with the magnetic axis sees the beacons separated 
by the unipolar induction charging period, instead of the pulsar rotation period. 

Here, we have reached the pulsating picture of the magnetosphere through a steady state 
approach by viewing the unipolar induction inside the pulsar and the magnetosphere outside 
as a coupled system. As a matter of fact, the steady state approach may not be adequate to 
describe a charge separated magnetosphere. We call to the attention that a charge separated 
plasma is a very unstable plasma due to the fact that any plasma flow (drift) is a current 
flow which generates magnetic field that drives drifts. For this reason, we might have to let 
go the lighthouse paradigm, and consider the time dependent model with Maxwell equations 


26 


coupled to the pi 

asma 

2013; 

Urpin 

I2OI2, 

to 

0 


asma equations flKomissarov 


2006 


Spitkovskvl I2OO6I: iTchekhovskov et al. 
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magnetosphere could be the cause of periodic beaming of electromagnetic waves. 
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10. Conclusions 

We have taken the view of interpreting the magnetic held B of the magnetosphere 
as stationary in space, and the electric held E as established by uploading the potential 
distribution of the pulsar surface to the magnetic held lines. The plasma velocity v is, 
therefore, the drift response in Eq. [2] which rotates with an angular velocity of the pulsar. 
The presence of a magnetosphere allows the pulsar surface charges be channelled to the 
held lines, and keeps the counter electric held in the pulsar interior low. This warrants 
the continuous operation of the unipolar induction that constantly pumps energy to the 
magnetosphere. In this scenario, the energetics of the magnetosphere is coupled to the 
unipolar induction of the neutron star. 

By introducing the gradient B and curvature plasma drifts and the presence of a toroidal 
held, we have removed the inconsistency between the space charge density and the curvature 
of the dipole-like closed ion held lines. Also, the inconsistency between introducing the far 
zone interstellar plasma boating potential to the near zone open held lines to generate an emf 
driven plasma how along open held lines and the E\\ = 0 plasma condition can be removed 
as well by drift velocities. Although plasma inertia and pressure are neglected in the force- 
free pulsar equation, they enter the magnetospheric dynamics implicitly through the E x B, 
gradient B, and curvature plasma drifts, which are the respective guiding center drifts of the 
basic gyro motion. Since the gradient B and curvature drifts are mass dependent, the closed 
region should be bounded at a distance within LC. 

We have devised a method to solve the pulsar equation analytically for the standard 
split monopole magnetosphere. The solutions are given in terms of three parameters, the 
separation function parameter /(T) = the axial current parameter A(T) = —and 
the separation constant m^. 

Recognizing that the magnetosphere is described by the helds of Eq. [T] under force-free 
approximation, and by the plasma drift velocities of Eq. |2] where plasma inertia and plasma 
pressure are prime factors, and these two descriptions are connected by the space charge 
densities of Eq. [101 we have presented a self-consistent description of the pulsar magneto¬ 
sphere. Through the coupled system between the unipolar induction inside the neutron star 
and the magnetosphere outside, and considering the unipolar induction rate exceeds the 
open held line Poynting hux, we have constructed a dynamic pulsating magnetosphere. In 
this pulsating model, the accumulated magnetic and plasma energies of the closed held lines 


can be released periodically by switching open the held lines either as plasma pressnre gets 
high enongh or as the toroidal magnetic held gets large enongh. As an alternative to the 
lighthouse paradigm, this pulsating magnetosphere ohers a diherent mechanism to generate 
pulsar beacons. 
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Fig. 2.— With = 5 corresponding to a| = 8 and = 3, the axial function Za(‘^) with 
exponentially decaying prohle is plotted as a function of the normalized axial distance 







Fig. 5.— With = 5 corresponding to a| = 6 and = 1, the axial function Za(<^) with 
exponentially decaying prohle is plotted as a function of the normalized axial distance 



for the equatorial region ion open fi 
015, 0.020, 0.025 with the same ordering 






Fig. 7.— With oo = 1 and Oi = —1.5, and = —17 and = 7, plus the condition 
of al = 0, the radial solution R{C) with a root within LC is plotted as a function of the 
normalized radial distance 7. 
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Fig. 9.— A set of poloidal flux function for the equatorial region ion closed field lines 

is illustrated with contour values C = 0.03, 0.04, 0.05 with the same ordering. These field 
lines extend to an intermediate range of ^ to assure < 0 as required by the space charge 
density. 









